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+^ Abstract 

O In this work, the area and edges spin wave calculations were carried out using the Heisenberg Hamilto- 

^—i nian and the tridiagonal method for the 2D ferromagnetic square lattice stripe, where the SW modes are 

characterized by a ID in-plane wave vector q^-. The results show a general and an unexpected feature that 
the area and edge spin waves only exist as optic modes. This behavior is also seen in 2D Heisenberg anti- 



ferromagnetic square lattice. This absence of the acoustic modes in the 2D square lattice is explained in |[T1 
^ by the fact that the geometry constrains for NN exchange inside the square lattice allow only optical modes. 

We suggest that this unexpected behavior of spin waves in the 2D square lattice may be useful in realizing 
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I. INTRODUCTION 



Spin waves in 2D magnetic systems are very interesting both experimentally [|2]-l6l and theoret- 
ically [TTl-fTSl. For example, these systems are relevant to our understanding of high temperature 
superconductors [|T4] - l20l . and are the basis of many technological applications of ultrathin ferro- 
magnetic films (e.g., magnetic memory and storage devices, switches, giant magnetoresistance, 
etc), as well as in the new promising field of spintronics (see [|2TI - l26l '). 

Many theoretical techniques has been used to study spin waves (SWs) in 2D and 3D Heisenberg 
magnets [10]. Some examples are the Holstein-PrimakofF (HP) method [1], the "boson mean- 
field theory" [|71 [H where Schwinger bosons are used to represent the spin operators, and the 
"modified spin-wave theory" [fTTl where the Dyson-Maleev transformation is used to represent the 
spin operators. Additionally, the semi-classical approaches [|27ll28l are widely employed for SWs 
at long wavelengths (or small wavevectors). In recent years, ultra-thin magnetic nano structures 
have been fabricated and studied extensively lE^ 1^91 - 1^ , for their SW dynamics. It is often useful 
to distinguish between the propagating SW modes and the localized SW modes. In a film geometry 
these are usually referred to as "volume" (or "bulk") modes and "surface" modes respectively. 
Thus, new theoretical studies are needed where the surfaces (and/or interfaces) are important and 
where the localized SW modes are considered in detail. 

In this work our aim is to study SW modes at low temperatures (compared to the Curie temper- 
ature Tc). We do this for ultra-thin ferromagnets with typically one atom thickness, specifically a 
2D stripe with a finite number of atomic rows (a nanoribbon). We employ the Heisenberg model 
assume, for simplicity, a square lattice. We study both the "area" SWs (that propagate across 
the stripe width) and localized "edge" SWs. This work is very interesting for its expected novel 
fundamental physics and promising application in magnetic devices. 

An operator equation-of-motion technique known as the tridiagonal matrix method Il33] - [36ll 
will be conveniently employed here to calculate the SW spectra and, in particular, to distinguish 
between the edge modes and area modes of the ferromagnetic square lattice nanoribbons. 

II. THEORETICAL MODEL 

The system initially under study is a 2D Heisenberg ferromagnetic stripe (or nanoribbon) in 
the xy-plane. We assume a square lattice with lattice constant a and we take the average spin 
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alignment of the magnetic sites to be in the z direction, which is also the direction of the applied 
magnetic field. The nanoribbon is of finite width in the y direction with atomic rows (labeled 
as n = 1, • • • , A^) and it is infinite in the x direction (-00 <^ 00). The position vector for each site 
is given by r = a{m,n,0), where m is an integer from -00 to 00, and n is the row number with 
n = 1,2, ■• • ,N (see Figure[T]). 



n=N 




FIG. 1. Geometry of a 2D Heisenberg ferromagnetic square lattice nanoribbon. The spins are in the xy- 
plane and the average spin alignment is in z direction. The nanoribbon is finite in y direction with atomic 
rows (n - \, - ■ ■ ,N). 

The total Hamiltonian of the system is given by 

the first term is the Heisenberg NN exchange term, the second term is the Zeeman energy term 
due to an applied field Hq, and the third term represents the uniaxial anisotropy. The summations 
over i and j run over all the sites, where the NN exchange Jij has the "bulk" value J when either 
i and j are in the interior of the nanoribbon, and a modified value when / and j are both at the 
edge of the nanoribbon (i.e., in row n = I or n = N). Similarly, for the uniaxial anisotropy term, 
we assume a value D when the site i is inside the nanoribbon, and a modified value for sites at 
the edge of the nanoribbon. 

To calculate the SWs for this system at low temperatures T <^ Tc where the spins are well 
aligned such that the thermal average of - S for each spin, we use the Holstein-Primakofl" (HP) 
transformation and follow similar procedures to UJ to express the total Hamiltonian in terms of 
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boson operators. We arrive to the expression 

^Total = Eq + Hs (2) 
where the constant term Eq is the energy of the ground state for the ferromagnetic system given by 



-^y.J'J-T.Di -8MbHoJ]s, (3) 



2^ ^ 



and the operator term has the following quadratic form 



where bj and bj are the creation and the annihilation boson operators. 

In order to diagonalize and obtain the SW frequencies, we may consider the time evolution 
of the creation and the annihilation operators b^j and bj, as calculated in the Heisenberg picture in 
quantum mechanics. In this case, the equation of motion [iJTl - HTI (using units with fi = I) for the 
annihilation operator bj is 

dbj(t) r 1 

= (guM, + (2S -l)Dj)bj(t)-^sY, JiJ (*.(') - */')) (5) 



where the commutation relation between b] and bj in 

[bi, b^j\ = 5ij, [b], b] = -6ij, [bi, b] = [b], b]] = 0. (6) 

was used, as well as the operator identity [AB, C] = A[B, C] + [A, C]B. 

The equation of motion for the creation operator b^j is easily obtained by taking the Hermitian 
conjugation of Equation ([5]), giving 

db](t) 1 

i-^ = -(gfiBHo + {2S - \)Dj)b]{t) + Z '''-i Vi^^^ - ^j^^n ■ 

The dispersion relations of the SWs (i.e., the energy or frequency versus wavevector) can now 
be obtained by solving the above operator equations of motion. The SW energy is related to the SW 
frequency using E = fio), and a Fourier transform for the operators from the time representation to 
the frequency representation is made: 

bj(x, a))e~''^'doj, 

CO 

Xoo 
b]{x, co)e-''^'dco. (8) 
CO 
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On substituting Equation ([8]) in Equations [5] and |7} we get 

( 1 "1 1 

cj - (embHo + (25 - l)Dj) -:zSj] J.j bj{co) + J^^^) = 0, 

( 1 "I 1 

oj + (g^iBHo + (25 - l)Dj) + E -^'-J ^)^'^^ ■^'.i^J(^) = 0- (9) 

V i J i 

Since the nanoribbon extends to +oo in the x direction, we may introduce a ID Fourier trans- 
form to wavevector along the x direction for the boson operators Z^j and bj as follows: 

bj(x,oj) = ybn{q,,oj)e''^^"'\ 

b]{x, a;) = 4= E ^nM'., co)e'<^""\ (10) 

where A^o is the (macroscopically large) number of spin sites in any row. The transformed operators 
obey the following commutation relations: 



\bn{qx,oS),bl{q^,io)\ = 5^^^'^. 



(11) 



By substituting Equation pO] ) into Equation (|9]) and rewriting the summations, we get the following 
set of coupled equations: 

L - {g^iBH^ + (25 - 1)D,) - ^5(27, + 7) + ^5 7,r(^,) j Z7w(^„ a;) 



+ JbN-\{qx,co) = fox n = N 

[co - {gfiBHo + (25 - 1)D) - ^5(47) + ^57r(^,)jz7„(^,,a;) 

+ Jibn+iiqx,oj) + bn-iiqx,oj)) = 
L - (siubHo + (25 - 1)D,0 - ^5(27y + 7) + JMqx)\biiqx,co) 



(12) 



for > n > 1 



+ ^SJb2(qx,co) = 



for n = I. 



The first and the third equations refer to the edges n = N and n = 1 for the nanoribbon system, 
and we have defined yiqx) = 2 cos(^v6')- Similar results can be found for the equations involving 
the creation operator b^j. 

The above coupled equations can conveniently be written in matrix form as 



(-QI + A)b = 0, 
(Q/ + A)b'^ = 0, 



(13) 
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where b and are x 1 column matrices whose elements are the boson operators b„(qx, co) and 
bl(q^, co). Also Q. = a>/S J is a dimensionless frequency. The second equation is redundant in that 
it does not give rise to any new physical modes, so we will therefore ignore it. Here / is the N xN 
identity matrix and A is the following tridiagonal N x N matrix: 



A = 



a, -I • 
-1 a -1 • 
-1 a -1 • 



•■■ a -1 
••• -1 a -1 
•■■ -1 a, ^ 
The following dimensionless quantities have been defined: 

_ 2(gfiBHo + (25 - 1)D,) - 5(27, + J) + S JMqx) 

2{giUBHo + (25 - l)D) - 5(47) + 57y(^,) 
' = 57 • 

It is convenient to denote a new matrix by A' = -Q.I + A, so 



(14) 



(15) 



where 



A' = 



a; -10 
-1 a' -1 
-1 a' -1 



a' -1 
-1 a' -1 



-1 a: 



a. 



-oj + 2(g//g//o + (25 - l)De) -S(2J, + 7) + 57,y(g,) 

57 

-io + 2igiUBHo + (25 - 1)D) - 5(47) + 57y(g,) 
57 



(16) 



(17) 



The new tridiagonal matrix A' may be separated into two terms, following the approach in [|35l[36l 
, as 

A' =Ao + A, (18) 



where 



a' -1 
-1 a' -1 
-1 a' -1 



a' -1 
-1 a' -1 
-1 a' 



(19) 



A = 



A 





(20) 



••■0 
•••0 
•••0 A ^ 

and the element A = a\-a' . In this way all the edge properties have been separated into the matrix 
A. The inverse of a finite-dimensional tridiagonal matrix with constant diagonal elements such as 
Ao is well known JSH |36l and can be expressed as 



(21) 



(1 -x2^+2)(^__^-l) 

Here x is a complex variable defined such that \x\ < 1 and x + x'^ = a'. 

On noting that A' = (Aq + A) = Ao(l + A~^\), the dispersion relations are obtained by the 
condition DM 33 that detA' = 0, which implies 



det(/ + Ao^A) = 



(22) 



Using the previous equations, the matrix M = (/ + A^'A) can next be written in a partitioned 
form [113: 



M = 










M2,l 








I 








Mn-i,n 


Mn,i 





Mn,n , 



(23) 
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where the nonzero elements of M can be written as 

Mij = 6ij + 6ijiA-\jA + 6nj{A-\n^, (24) 



and the determinant of M, which is required for Equation ^22^ , can be calculated to give 

det(M) = (Mu)'-(Mi,;v)' 

-A 



After some more algebraic steps the condition for det(A/) = can be written as 

[(l - x^^^^) (x - X-') + (x^ + x^^ - x^^^^ - l) a] 

+77(2x^"^-;c^-1-x^^^)A = 0, (26) 

where 77 = +1. This result is formally similar to an expression obtained in the study of finite 
thickness ferromagnetic slabs [|34l . Also, it follows by analogy with previous work [|34l [35l that 
the solutions with \x\ = 1 correspond to the area modes (those propagating across the width of the 
stripe) while those with \x\ < 1 correspond to the localized edge modes. 

Before presenting a general analysis of Equation ( |26l ) for the SW frequencies in Section III we 
next examine some special cases. 

A. Special case of A'^ ^ 00 

It is of interest to study the behavior of the model in the special case when the ribbon (or stripe) 



becomes very wide. This is the limit of ^ 00 in Equation ( [26| ). Since 1^1 < 1 for edge modes, 
the terms of order ^ as A^^ ^ 00, giving for these modes 

(l+xA) = ^x = -^. (27) 

This is the same formal expression as obtained in the case of a semi-infinite Heisenberg ferromag- 
net [|3ll35l when N ^ cx>. 



For the edge mode localization condition \x\ < 1 to be satisfied in this special case. Equation 27 

implies |A| > 1, which gives two possibilities. The first one is A > 1 and the second one is A < -1. 

From the definitions of a'^ and a' we have 

^ 2(25 - 1)(D, - D) - 5(27, - 37) + 5(7, - 7)y(g,-) 

57 
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which depends on the physical parameters D^., D, S , Je and 7 of the ferromagnetic stripe and on 
the wavevector component q^. Since y(<?x) = 2 cosiq^d) has its maximum when cosiq^a) = 1 and 
its minimum when cos(^v^*) = -1, it follows that the minimum value A^m and the maximum value 
^max for A correspond to 

_ 2(25 - 1)(D, -D)-S (27, - 37) + 25 (7, - 7) 

In simple cases we might have Dg = D i.e., the edge perturbation to the anisotropy is negligible. 
Then, denoting the ratio between the edge exchange and area exchange by r = 7^/7 > 0, we have 

The two cases A > 1 and A < - 1 give the following ranges for r. 

5 6 
r<- r>-. (28) 

These are the ranges of the ratio exchange r = J^/J for which the edge modes exist at some q^ 
value in this special case (A^ oo) for the square lattice. 

The conditions are modified if edge perturbations in the anisotropy are included. 



B. Case of large finite A'^ 

Another interesting case is when becomes sufficiently large that the two solution and x' 



for 77 = +1 of Equations ( [26| ) can be obtained by an iterative approach used in reference [|34l . Since 
\x\ < 1 for edge modes, all terms of order x^ in Equations ([26]) are small for sufficiently large N 



and the two solution and x become closer to the solution xq = -A ^ for the special case of 
A^ ^ 00. 



To use a first order iteration for Equations ([26]) we must rewrite them in the forms x- = F-{xo). 
We get the following two first order iteration approximate solutions 

-1 



■^ap/a 

where Xap is equal to 



(29) 



and Xam is equal to 



;cf ^3 + Axf - ^> - Axf - Ax^^' + 2Ax^^' - Ax^-' + A. 



0.96 

0.94 

0.92 

0.90- 

0.88- 

0.86- 

0.84- 

0.82- 

0.80- 

0.78- 

0.76- 



A = -1.05 



— I— 

20 



A = -1.1 



A = -1.3 



— 1 ' 1 ' r- 

40 60 80 

N number of rows 



100 



FIG. 2. Calculated values of and x for several values of A < -1 and for N from 10 to 100. 



This approximation is valid provided |A| > 1. Figure [2] shows calculated values of and x' using 
illustrative values of A chosen as -1.05, -1.1 and -1.3 when has values from 10 to 100. The 
iterative approach is found to work well for N > 10 in this case. 



ni. NUMERICAL CALCULATIONS 



More generally, the dispersion relations can be obtained by solving Equation ( [26] ) using a nu- 
merical calculation for any finite A^. The number of rows and the value of A are first substituted 



in the Equations ( [26| ), and then the polynomial equations are solved for x which can be used to 
obtain the dispersion relations. Since the solutions for x may have complex roots, one of the ways 
to solve such equations is to use Laguerre's method for finding the roots [|44ll45l of polynomials. 
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By rearranging Equations ( |26l ) we have: 



First polynomial 

-Ax^ -x^ + Ax+l=0 (30) 
Second polynomial 

^2N.4 ^ ^^2N.3 _ ^2N.2 _ ^^2N.i ^ _ 2^^N.2 ^ ^^N 

-Ax^ -x^ + Ax+l=0 (31) 

Both polynomials are of degree IN + 4, and they are applicable for all equal to 3 and above. 
We note that there is a special case when N = 3, since the two power indices 2N + 1 and A'^ + 4 
become equal to 7. The obtained values for x must satisfy the conditions for physical SW modes 
mentioned earlier. The edge SW modes are localized on the edge and their amplitudes decay 
exponentially inside the nanoribbon. This requires that x must be real and less than 1 for edge 
modes. The area modes are oscillating waves inside the nanoribbon, and so 

X e R and \x\ < 1 for edge modes 

X = e'*"' and \x\ < 1 for area modes (32) 

In the previous section the ranges for r for edge modes to exist were obtained algebraically for 
the special cases N ^ oo and large A^. We can use numerical calculations to obtain the ranges of 
A for smaller that satisfies the conditions for which both edge modes and area modes exist. 

For that purpose, a Fortran program was written to solve the two polynomials using Laguerre's 
method where two subroutines zroots and laguer are adapted from fl44ll . 

The values for minimum positive (P) and maximum negative (N) of A for even (E) and odd 
(O) rows number N that satisfy ( [32] ) are computed from the first polynomial (F)|30]and from the 
second polynomial (S)[3T]and displayed in figure[3]for edge modes and in figure|4]for area modes. 

Figure |3] shows the behavior of minimum positive for A that satisfies the conditions ( [321 ) for the 



existence of edge modes. It is clear from the figure, that in the range of rows number A^ < 20, 
the minimum positive of both odd rows number of first polynomial [30] and even rows number 
for second polynomial [3T] are approximately the same and are exponentially decaying to a nearly 
constant value of 1.02. In the same rows number range, the minimum positive of even rows for 
first polynomial [30] and of odd rows for second polynomial [3T]is nearly constant and equal to 0.95. 
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FIG. 3. The values for minimum positive (P) and maximum negative (N) of A for even (E) and odd (O) 



rows number N that satisfies edge modes ( 32 1, are computed from the first polynomial (F) 30 and from the 
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FIG. 4. The values for minimum positive (P) and maximum negative (N) of A for even (E) and odd (O) 
rows number N that satisfies area modes ( |32l ) are computed from the first polynomial (F) 30 and from the 



second polynomial (S) 31 
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After equal to 20, the minimum positive of both polynomials is independent of stripe (ribbon) 
width, (i.e. rows number), and it is also independent of rows number parity, whether even or odd. 
As the rows number increases, the minimum positive is convergent to an approximate constant 
value of 1.02. 

The same Figure |3] displays the maximum negative of A for first polynomial 30 and second 
polynomial |3T] Here, the maximum negative shows different behavior from that of the above 
minimum positive as it is independent of rows number parity for both polynomial. In the range of 
rows number A'^ < 20, as in case of minimum positive, the maximum negative of first polynomial 



30 is nearly constant and equal to -0.96, while the maximum negative of second polynomial 31 



exponentially increases to a nearly constant value of 1.02. As in the case of the minimum positive, 
after equal to 20, the maximum negative of both polynomials is independent of stripe (ribbon) 
width (i.e. rows number), and rows number parity, whether even or odd. As increases, the 
maximum negative value is convergent to approximately constant value equal to -1.02. 

The conclusion from Figure[3]is that edge modes, in small rows number A'^ < 20, are dependent 
on both the stripe width and the rows number parity. This is an indication for the interaction 
between the two edges in the small range of rows number. As A^ increases above 20, the edge 
modes become independent on both the stripe width and the rows number parity. In this case both 
minimum positive and maximum negative of both polynomials are stripe width and rows number 
parity independent. This is an indication for disappearing of the interaction between the two edges 
after A^^ = 20. That behavior agrees with result for the special case of A^^ become large, as discussed 
above. Also, it is noted in the range of A^ larger than 20, that the difference between minimum 
positive and maximum negative is nearly constant and independent of the stripe width and rows 
number parity. 

Figure|4]shows the behavior of minimum positive and maximum negative for A that satisfies the 



conditions ( [32] ) for area modes. It is clear that both minimum positive and maximum negative are 
independent of rows number parity. While they depend on the stripe width, in the range of rows 
number from 5 to 10 the two values are constant. After A'^ = 10, the value of minimum positive 
is increasing linearly with the stripe width, wliile the value of maximum negative is decreasing 
hnearly with the stripe width. The difference between minimum positive and maximum negative 
is increasing as the stripe width increases and it is mostly independent on the rows number parity. 

To obtain the dispersion relations for the above system, a Fortran program was written to solv- 
ing the first polynomial [30] and the second polynomial [3T]by Laguerre's method using two subrou- 
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tines zroots and laguer adopted from (44], which are used before. The values of physical pa- 
rameters for calculating these dispersion relations are chosen as follow: 5 = l,7=l,Z) = Dg = 
and giUsHo = 0.3 J. The chosen value for the ratio between the edge exchange and area exchange is 



equal to r = Jg/J = 0.04, which is satisfies the existence condition 28 for edge mode. The chosen 
values of q^a are run from to ;r corresponding to the first Brillouin zone center and boundary 
respectively. 



IV. SW DISPERSION RELATIONS 

The numerical results for calculating the dispersion relations of 2D square lattice using the 
above algorithm and physical parameters are displayed in Figures [5]j8] These figures show plots 
of SW frequency in terms of the dimensionless quantity co/S J as a function of the dimensionless 



wavevector q^a for various numbers of rows. The two polynomials 30 and 31 have even power in 
X, and therefore the area and edge modes are symmetric about co/S J = 0. As a result, we have 
chosen to show only the positive frequency branches. 

All the figures display similar general features for ferromagnetic 2D stripes. 

Since every polynomial gives area modes, the number of area modes are equal to twice the 
number of rows A^. These area modes are upper bounded by the frequency obtained when the value 
of X is equal to = 1, and lower bounded by the frequency obtained when the value of x is equal to 
= -1, as shown in the Figures. In the upper of area modes, some spin wave frequencies cross each 
other. As A'^ increases, the number of areas modes inside their boundary increases twofold, which 
then merge into an areas modes continuum. 

The figures show that in all cases that there are two optic edge modes appearing above the area 
modes region, these two edge modes look like extension for their counterpart area modes. As A'^ 
increases, the diff"erence between the two edge modes is decreasing which is seen too for their 
counterpart area modes. 



V. DISCUSSION 

The dispersions relations of area and edges spin waves, and the eff'ect of the stripe width on 
them for 2D ferromagnetic square lattice stripes, have been studied using the tridiagonal method. 
The result shows the same unexpected feature: the area and edge spin waves only exist in optic 
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Dispersion relation for N = 3 




FIG. 5. Area and edge spin waves modes (in units of S J) plotted against the wavevector q^a for stripe with 
width N = 3, where x = I and x = -I are the upper and lower boundary for area modes. 
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FIG. 6. Area and edge spin waves modes (in units of S J) plotted against the wavevector q^a for stripe with 
width N = 4, where x = I and x = -I are the upper and lower boundary for area modes. 
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FIG. 7. Area and edge spin waves modes (in units of SJ) plotted against the wavevector q^a for stripe with 
width N = 1, where x - \ and x - -\ are the upper and lower boundary for area modes. 




FIG. 8. Area and edge spin waves modes (in units of S J) plotted against the wavevector q^a for stripe with 
width N = where ;c = 1 and ;c = -1 are the upper and lower boundary for area modes. 
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modes. This behavior is also seen in 2D Heisenberg antiferromagnetic square lattice experimental 
and theoretically [|46] - l48l . the absence of the acoustic modes could be explained that the square 
lattice support only optics modes, which need more studies to be completely understood. 

Our conclusion, that the unexpected behavior of spin waves in the 2D square lattice of existence 
in only optic modes if included in the HTS theories may lead to an explanation for HTS. Since it 
is known that HTS is linked to 2D square antiferromagnetic lattice, and we expect that the Optic 
spin wave could mediate the electrons using their spin degree of freedom and Pauli exclusion 
principle for the formation of cooper pair with much less energy than cooper pair created by 
phonon mediated electrons using their electric charge [[T5ll49] - l53]| . 
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